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ABSTRACT 


Analysis  of  laminar,  transitional,  and  turbulent  boundary  layers  on 
a  sharp  cone  at  zero  angle  of  attack  under  hypersonic  perfect  gas  condi¬ 
tions  is  presented.  The  governing  boundary -layer  equations  were 
numerically  integrated  on  a  digital  computer  using  a  marching,  iterative, 
implicit,  finite -difference  method.  The  turbulent  boundary  layer  is  ana¬ 
lyzed  using  a  two-layer  eddy  viscosity  model  with  a  constant  turbulent 
Prandtl  number,  and  the  transition  region  is  treated  through  an  eddy 
viscosity-intermittency  factor  approach.  Comparison  with  experi¬ 
mental  data  reveals  that  the  present  approach  yields  an  accurate  de¬ 
scription  of  laminar,  transitional,  and  fully  turbulent  heat  transfer  on 
sharp  cones  in  the  Mach  number  range  from  5  to  10  under  cold  wall  con¬ 
ditions.  Verification  of  the  detailed  calculation  of  transitional  and  turbu¬ 
lent  boundary -layer  structure  under  hypersonic  conditions  must  await 
further  experimental  investigations.  It  appears  at  the  present  time  that 
the  eddy  viscosity  approach  to  calculation  of  transitional  and  turbulent 
boundary  layers  may  indeed  be  applicable  and  accurate  under  hypersonic 
conditions. 
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SECTION  I 
INTRODUCTION 


The  successful  design  and  operation  of  a  reusable  lifting  entry  space¬ 
craft  depends  partly  upon  the  design  of  the  thermal  protection  system.  In 
turn,  the  thermal  protection  system  design  depends  primarily  on  the 
thermodynamic  environment.  Peak  heat-transfer  rates  and  duration  of 
the  heating  are  the  most  influential  thermodynamic  parameters.  In 
general,  the  peak  heat -transfer  rate  determines  the  thermal  protection 
system  materials,  and  the  duration  of  heating  and  the  integrated  energy 
available  determines  the  mass  of  the  thermal  protection  system. 

Boundary-layer  transition  to  turbulence  is  well  known  to  have  a 
major  influence  on  the  magnitude  of  lifting  entry  body  heating  rates. 

Lifting  entry  trajectories,  characterized  by  a  gradual  variation  in  altitude 
and  velocity,  magnify  the  influence  of  the  transitional  zone  on  local  heating 
histories.  The  gradual  variation  of  local  conditions  induces  the  location 
for  transition  onset  to  move  forward  at  a  slow  rate.  At  a  fixed  location 
on  the  body,  the  fully  turbulent  heat -transfer  rate  decreases  with  time, 
causing  peak  turbulent  heating  to  occur  near  the  inception  of  fully  turbu¬ 
lent  flow.  With  different  assumptions  of  the  development  time  necessary 
to  achieve  fully  turbulent  flow,  significant  variances  in  local  heating 
histories  are  possible.  Hence,  accurate  predictions  of  peak  heating  rates 
and  total  integrated  heating  loads  require  a  realistic  model  of  the  transi¬ 
tion  zone.  The  assumption  of  fully  turbulent  flow  at  a  point  where  the  flow 
is  actually  transitional  can  result  in  an  overweight  design  of  the  thermal 
protection  system,  which  cannot  be  tolerated  because  of  the  low  ratio  of 
payload  to  gross  weight  characteristic  of  lifting  entry  vehicles.  This 
subject  has  been  the  topic  of  a  recent  paper  by  Masaki  and  Yakura 
(Ref.  1). 

Prediction  techniques  for  the  nonsimilar  fully  turbulent  hypersonic 
boundary  layer  utilizing  large  digital  computer  codes  have  become  avail¬ 
able  during  the  past  few  years,  e.  g. ,  Smith  and  Cebeci  (Ref.  2), 

Patankar  and  Spalding  (Ref.  3),  Bushnell  and  Beckwith  (Ref.  4),  Herring 
and  Mellor  (Ref.  5),  Martellucci,  Rie,  and  Sontowski  (Ref.  6),  and 
Mayne  and  Dyer  (Ref.  7).  Comparison  of  results  from  the  above  analyses 
with  experimental  data  have  shown  good  agreement  for  fully  turbulent 
hypersonic  (up  to  M,,,*  10)  boundary  layers  on  sharp  and  slightly  blunted 
cones  as  well  as  sharp  flat  plates.  These  results  suggest  that  the  eddy 
viscosity  model  employing  the  Prandtl  mixing  length  is-  indeed  adequate 
for  prediction  of  hypersonic  turbulent  boundary -layer  flows  over  simple 
bodies  under  cold  wall  conditions  and  may  be  used  with  some  confidence. 
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With  regard  to  prediction  of  transitional  boundary -layer  flows, 
little  has  been  accomplished.  Most  of  the  references  above  employ  an 
instantaneous  step  transition  from  fully  laminar  to  fully  turbulent  flow 
at  either  a  specified  body  location  or  local  Reynolds  number  along  the 
boundary -layer  edge.  An  illustration  of  this  type  of  treatment  is  given 
in  the  sketch  below. 


Local 

Heat- 

Transfer 

Rate 


It  is  readily  apparent  that  the  instantaneous  transition  approach  is  not  a 
valid  model  for  the  transition  region.  In  order  to  provide  a  suitable 
progressive  transition  model,  Masaki  and  Yakura  (Ref.  1)  propose  an 
intermittency  factor  phenomenological  description  of  the  transition 
region  based  on  the  experimentally  observed  fact  that  transition  in  a 
boundary  layer  is  characterized  by  the  intermittent  appearance  of 
turbulent  spots  which  grow  and  move  downstream  with  the  fluid. 
Although  this  approach  appears  to  have  merit,  Masaki  and  Yakura  pre¬ 
sent  no  results  for  the  progressive  transition  model  relative  to  experi¬ 
mental  data.  They  state  that  the  primary  purpose  of  their  paper  is  to 
encourage  further  studies  on  this  subject. 
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The  present  work  represents  an  adaptation  of  Masaki  and  Yukura's 
intermittency  factor  description  of  the  transition  region  coupled  with 
implicit  finite -difference  numerical  solution  of  the  governing  turbulent 
boundary -layer  equations  for  the  case  of  hypersonic  flow  over  a  sharp 
cone  at  zero  angle  of  attack.  Since  the  current  state-of-the-art  in  turbu¬ 
lent  boundary -layer  analysis  under  hypersonic  conditions  is  very  un¬ 
certain,  comparison  of  theoretical  predictions  with  experimental  data  is 
presented  in  order  to  establish  the  basic  validity  of  the  present  approach. 
In  general,  the  intermittency  factor  treatment  used  herein  appears  to 
yield  a  reasonably  accurate  description  of  transitional  he'ating  on  sharp 
cones  in  hypersonic  flow  based  on  comparisons  with  experimental  data. 
This  finding  says  very  little,  however,  about  the  detailed  prediction  of 
transitional  and  turbulent  boundary -layer  structure  under  hypersonic 
conditions  which  must  await  further  experimental  investigation. 

The  favorable  results  of  the  present  study  indicate  that  the  inter¬ 
mittency  factor  approach  coupled  with  numerical  solution  of  the  turbu¬ 
lent  boundary -layer  equations  may  be  a  powerful  tool  in  analysis  tech¬ 
niques  for  transitional  heating  on  lifting  body  configurations.  Much 
work  remains  to  be  done  in  this  area. 


SECTION  II 

THEORETICAL  MODEL 

The  theory  used  herein  is  based  upon  numerical  solution  of  the  non¬ 
similar,  compressible  boundary -layer  equations  written  in  transformed 
variables  using  a  marching,  iterative,  implicit  finite -difference  pro¬ 
cedure.  The  turbulent  transport  terms  in  the  boundary -layer  equations 
are  treated  using  a  two-layer  eddy  viscosity  model  based  on  Prandtl's 
mixing  length  in  conjunction  with  a  modification  of  van  Driest' s  analysis 
for  the  near-wall  region.  The  transition  zone  is  represented  by  an 
intermittency  factor  approach  coupled  to  the  eddy  viscosity  model. 


2.1  GOVERNING  BOUNDARY-LAYER  EQUATIONS 


The  present  analysis  employs  the  compressible  turbulent  boundary- 
layer  equations  in  terms  of  time -averaged  mean  flow  quantities.  With 
the  following  assumptions  for  the  Reynolds  stress  and  heat  transfer 

-  p"  ITF7  =  e  (1) 


-  p  v'n'  =  k 


3T  _  k  3tT 

3y  "  ^  3y 


(2) 
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and  the  definition  of  the  turbulent  Prandtl  number 

Cde 

Prt  =  -f"  «) 

the  governing  compressible  boundary -layer  equations  describing  the 
flow  about  two-dimensional  and  axisymmetric  bodies  at  high  Reynolds 
number  and  constant  pressure  across  the  boundary  layer  are  (Refs.  2, 
8,  9,  and  10): 


CONTINUITY 


3_ 

3X 


(p  V  rt,j) 


=  0 
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where 


(7) 


and  the  usual  expressions  for  the  mean  and  fluctuating  parts  of  the 
variables  are  used,  e.  g. , 


p  =  p  +  p'  (8) 

The  exponent  j  is  equal  to  zero  in  a  planar  case,  and  equal  to  unity  for 
an  axisymmetric  case.  The  intermittency  factor.  If,  and  the  eddy 
viscosity,  e  ,  are  discussed  in  greater  detail  later  in  this  section.  The 
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/ 

static  pressure  variation,  p(x),  is  regarded  as  input  to  the  problem  and 
hence  a  known  quantity  as  is  customary  in  boundary -layer  analyses. 

The  coordinates  are  a  curvilinear  system  in  which  x  is  distance  along 
the  surface  measured  from  the  stagnation  point  or  leading  edge.  The 
dimension  y  is  measured  normal  to  the  surface.  Within  the  boundary 
layer,  the  velocity  components  in  the  x-  and  y-direction  are  u  and  v, 
respectively.  The  body  radius  (for  an  axisymmetric  body)  is  r^.  See 
Fig.  1,  Appendix  I,  for  clarification  of  nomenclature  relative  to  sharp 
cone  geometry  which  is  of  sole  interest  in  the  present  work. 

If  subscript  w  denotes  wall  and  e  denotes  outer  edge,  the  associ¬ 
ated  boundary  conditions  on  the  above -defined  equations  are: 


MOMENTUM 

u(x,y=0)  =  0  (9) 

V(x,y=0)  =  0  (10) 

lim  u(x,y)  =  U  (x)  (11) 

y  -*  oo  e 

ENERGY 

T  (x,y=0)  =  Tw(x)  (12) 

lim  T  (x,y)  =  Te(x)  (13) 

y  ■*  « 


which  reflect  the  requirements  of  no  slip  or  mass  injection  (blowing  or 
suction)  at  the  wall  as  well  as  a  prescribed  wall  temperature  variation. 
The  outer  edge  velocity,  Ue,  and  static  temperature,  Te,  must  be  de¬ 
termined  from,  say,  an  inviscid  analysis  consistent  with  the  imposed 
static  pressure  distribution. 

In  order  to  dose  the  set  of  governing  equations,  it  is  further 
assumed  that  the  gas  is  thermally  and  calorically  perfect  air  having  a 
constant  specific  heat  ratio  y  =  1.  40  and  obeying  the  equation  of  state 

p  -  pRT  (14) 
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where  R  =  1716  ft^/sec^-  °R.  Hence,  under  this  assumption 

h  =  CpT  (15) 

where  Cp  =  6006  ft^/sec^  -  °R,  The  laminar  viscosity,  ^ ,  is  taken  to 
obey  Sutherland's  law 

£ _ _  W  +  198.6  /f  \  3/2 

yref  ~  T  +  198.6  l  TreT /  <16) 

where  T  must  have  units  of  °R  and  subscript  ref  denotes  a  reference 
condition.  The  laminar  Prandtl  number,  Pr,  is  taken  to  be  a  constant 
value  of  0.  71  across  the  entire  boundary  layer. 


2.2  TURBULENT  TRANSPORT  LAWS 

The  shear  stress  in  a  turbulent  boundary  layer  is  treated  herein  by 
the  use  of  a  two-layer,  inner-outer  model  using  Prandtl's  mixing -length 
hypothesis  and  a  modification  of  van  Driest' s  analysis  for  the  near-wall 
region.  This  results  in  a  continuous  distribution  of  the  shear  stress 
from  the  laminar  value  at  the  wall,  through  the  fully  turbulent  region, 
reaching  zero  at  the  outer  edge  of  the  boundary  layer.  The  energy 
transport  in  a  turbulent  boundary  layer  is  treated  in  this  work  through 
the  incorporation  of  the  eddy  conductivity,  *,  into  a  turbulent  Prandtl 
number,  Pr^,  according  to  Eq.  (3). 

Following  the  usual  treatment  of  the  boundary-layer  equations  for 
a  turbulent  two-dimensional  boundary  layer,  the  effective  shear  stress 
may  be  written  as 

T  =  v  —  -  pFv1'  (17) 

where  the  term  -  p  xFv'-  is  called  the  Reynolds  stress  and  represents 
the  shear  stress  introduced  by  the  turbulence.  Using  Prandtl's  mixing 
length  hypothesis  for  the  Reynolds  stress  results  in 

-  p  u'v'  =  p  l* 

where  J2*  is  called  the  mixing  length.  Comparison  of  Eqs.  (1)  and 
(18)  reveals 


3U 

ay 


(18) 


e 


-  2  3U 

P  i*  jjy 


(19) 
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which  is  the  definition  of  the  eddy  viscosity,  e  ,  for  use  in  the  govern¬ 
ing  boundary -layer  equations. 

After  Escudier  (Ref,  11),  Patankar  and  Spalding  (Ref.  3)  recommend 
the  following  variation  of  the  mixing  length,  i*,  across  the  boundary 
layer 

t*  =  k*.y,  for  0  <  y  *  Xy^/k*  (20) 

l*  =  Xy£,  for  xyA/k*<  y  (2i) 

where  the  values  for  the  various  numerical  constants  are  taken  to  be 
k,jc  =  0.435  and  X  =  0.  09.  The  value  of  y  at  the  point  where  the  velocity 
in  the  boundary  layer  is  equal  to  0.99  of  the  velocity  at  the  boundary- 
layer  outer  edge  is  used  to  define  the  distance  y The  above  choices 
follow  Patankar  and  Spalding  and  result  in  good  agreement  with  the 
mixing-length  model  proposed  by  Maise  and  McDonald  (Ref.  12)  for 
compressible  turbulent  boundary  layers. 

By  analogy  with  Stokes'  solution  for  an  infinite  flat  plate  undergoing 
simple  harmonic  motion  parallel  to  itself  in  an  infinite  fluid,  van  Driest 
(Ref.  13)  concluded  that  in  the  vicinity  of  a  wall  the  shear  stress  in  a 
turbulent  fluid  should  be  of  the  form 


t  = 


3U 


+  p  k*2y2 


f  _ _ \m 

2 

I  \7 

1  -  exp 

[Vtw  p 

li  a*  I 

ay 

\  IJ 

I  / 

(22) 


which  results  in  an  exponential  damping  of  the  turbulent  part  of  the 
shear  stress  as  the  wall  is  approached  and  yields  exactly  the  laminar 

shear  stress  form,  r  =  V  -r— ,  at  the  wall.  Although  Eq.  (22)  was  origi- 

oy 

nally  developed  for  incompressible  flow,  it  is  applied  in  the  present 
work  to  compressible  flow  by  application  of  the  suggestion  by  Patankar 
and  Spalding  (Ref.  3)  that  the  local  value  of  shear  stress  be  used  instead 
cf  the  wall  value  as  originally  recommended  by  van  Driest  (Ref.  13). 

The  constant  A*  is  taken  to  be  26.0  following  van  Driest  (Ref.  13). 


The  turbulent  Prandtl  number  as  defined  by  Eq.  (3),  i.  e. , 


LnE 

P*t  =  < 


is  physically  a  measure  of  the  ratio  of  eddy  viscosity  to  eddy  conduc¬ 
tivity,  that  is,  the  ratio  of  the  turbulent  transport  of  momentum  to  the 
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turbulent  transport  of  heat.  Since  the  flow  in  the  outer  region  of  a  turbu 
lent  boundary  layer  shows  some  similarity  to  a  turbulent  wake  flow,  one 
may  advance  arguments  that  a  realistic  formulation  of  turbulent  Prandtl 
number  requires  a  separate  expression  for  the  inner  and  outer  region 
just  as  in  the  eddy  viscosity  formulation.  A  recent  review  by  Rotta 
(Ref.  14)  indicates  that  the  turbulent  Prandtl  number  varies  in  magni¬ 
tude  from  approximately  0.  90  near  the  boundary -layer  outer  edge  to 
approximately  2. 0  near  the  wall.  Much  work  remains  to  be  done  in  de¬ 
fining  the  turbulent  Prandtl  number  distribution  as  the  review  by  Rotta 
points  out.  Hence,  in  the  present  work  the  turbulent  Prandtl  number  is 
taken  to  remain  constant  at  0.  90  across  the  entire  boundary  layer  as 
recommended  by  Patankar  and  Spalding  (Ref.  3).  This  choice  of  turbu¬ 
lent  Prandtl  number  is  consistent  with  the  analyses  of  Patankar  and 
Spalding  (Ref.  3),  Martellucci,  Rie,  and  Sontowski  (Ref.  6),  and  Mayne 
and  Dyer  (Ref.  7);  Smith  and  Cebeci  (Ref.  2)  use  a  value  of  unity. 


The  constraint  used  to  define  the  end  of  the  inner  region  and  the 
beginning  of  the  outer  region  is  the  continuity  of  the  eddy  viscosity. 
From  Eq.  (22),  the  eddy  viscosity  in  the  inner  region  is 


*  P 


k,y 


!’• 


exp 


2 


(23) 


whereas  for  the  outer  region 


-22 
pa2  y, 


3_U 

ay 


(24) 


with  the  constants  k*.  A*,  A,  and  y ^  as  defined  previously.  Note  that 
the  damping  term  in  Eq.  (23)  reflects  the  application  of  the  local  shear 
stress  as  opposed  to  the  wall  shear  stress  of  Eq.  (22)  as  discussed  pre 
viously.  From  the  wall  outward,  the  expression  for  the  inner  eddy  vis¬ 
cosity  applies  until 


ei  '  eo  (25) 

at  which  time  the  outer  eddy  viscosity  takes  over.  A  schematic  of  this 
variation  in  terms  of  the  mixing  lengths  is  shown  below. 
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2.3  TRANSITION  ZONE  DESCRIPTION 

The  distance  between  transition  onset  at  x  =  Xj-  and  the  beginning  of 
fully  turbulent  flow  farther  downstream  at  x  =  Xrp  is  denoted  herein  as 
the  transition  zone.  In  the  present  work,  which  is  restricted  solely  to 
analysis  of  sharp  cone  flows,  no  attempt  is  made  to  predict,  correlate, 
or  otherwise  a  priori  infer  the  location  of  transition  onset  at  X^.  In¬ 
stead,  the  experimentally  indicated  beginning  of  the  transition  region 
is  used  for  the  onset  point  X|.  based  on  the  particular  sharp  cone  heat- 
transfer  data  under  examination. 

Repeatability  of  gross  characteristics  of  the  transition  region  such 
as  a  nondimens ional  length  have  been  observed  in  experimental  investi¬ 
gations.  Potter  and  Whitfield  (Ref.  15)  show  that  the  transition  zone 
local  Reynolds  number  Ree>^x  (where  Ax  =  Xx  -  X^)  has  a  dependence 


9 


AEDC-TR-70-210 


on  both  transition  onset  local  Reynolds  number  Re6j  x^.  the  l°cal 

Mach  number  Me.  In  the  present  work,  the  extent  of  the  transition  zone 
will  be  represented  by  a  transition  zone  local  Reynolds  number  ratio 

Ree  XT^Ree  Xt  which  is  related  to  Ree,  Ax  hy 

Ree,XT  ,  ^ee,Ax 

1  +  Ree,Xt  <26> 

where  Ree  x-p  is  the  local  Reynolds  number  at  the  inception  of  fully 

turbulent  flow.  The  onset  Of  transition,  Xt,  in  the  present  work  is  deter¬ 
mined  from  examination  of  experimental  data  as  discussed  previously. 
The  end  of  transition  is  taken  herein  to  coincide  with  the  location  of 
experimentally  determined  maximum  heating.  A  schematic  of  this 
nomenclature  is  shown  below. 


Heat 

Transfer 


Surface  Distance  along  Body 
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Masaki  and  Yakura  (Ref.  1)  present  a  comprehensive  survey  of  un¬ 
classified  experimental  data  showing  transition  on  flat  plates  and  cones 
over  a  range  of  local  Mach  number  from  2  to  15.  Based  on  this  survey 
(which  is  given  in  Figs.  4  and  5  of  their  paper)  they  suggest  that  a  value 
for  Ree>  Xj/^ee,  Xj.  °f  approximately  2  appears  reasonable  at  high  super¬ 
sonic  speeds.  The  value  of  two  is  not  a  universal  number,  and  a  great 
deal  of  work  remains  to  be  done  before  this  boundary -layer  character¬ 
istic  is  fully  understood.  The  effects  of  wall  cooling  and  pressure  gradi¬ 
ent  are  moot  points  at  this  time.  It  should  be  noted,  however,  that 
Murphy  and  Rubesin  (Ref.  16)  found  that  Re6j  Xrp/Ree,  555  2  based  on 

flight  test  results  from  the  Mark -2  blunt-nosed  ballistic  reentry  vehicle. 
For  the  present  work 

R^.Xj  . 

Reo.Xt  "  (27) 

following  Masaki  and  Yakura  (Ref.  1).  As  more  hypersonic  data  become 
available,  this  factor  of  2  may  well  require  readjustment.  For  flows 
with  pressure  gradient  and  variable  wall  temperature  the  factor  of  2  may 
not  be  applicable.  All  that  can  be  said  for  its  present  use  is  that  it 
appears  reasonable  based  on  available  sharp  cone  and  flat  plate  data  from 
supersonic  and  hypersonic  wind  tunnel  flows. 


2.4  INTERMITTENCY  FACTOR 


In  order  to  provide  a  suitable  progressive  transition  model  as  opposed 
to  an  instantaneous  transition  model,  the  present  analysis  utilizes  the 
phenomenological  theory  of  the  transition  zone  developed  by  Emmons 
(Ref.  17)  which  is  based  upon  the  experimental  observation  that  transition 
in  a  boundary  layer  is  characterized  by  the  intermittent  appearance  of 
turbulent  spots  which  grow  and  move  downstream  with  the  fluid.  In  the 
transition  zone,  the  intermittent  appearance  of  turbulent  spots  and  their 
growth  as  they  move  downstream  produces  alterations  of  laminar  and 
turbulent  flow  with  a  gradual  increase  in  the  mean  turbulence  level.  At 
any  point  in  the  transition  zone,  the  probability  of  turbulent  flow  is  pre¬ 
scribed  by  a  model  based  on  the  residence  times  of  intermittent  turbulent 
spots.  Dhawan  and  Narasimha  (Ref.  18)  modified  Emmons'  original  formu¬ 
lation  of  the  probability  of  turbulent  flow  to  account  for  a  localized  region 
of  breakdown  about  x  =  Xf  in  conjunction  with  a  single  universal  intermit  - 
tency  distribution.  The  Dhawan  and  Narasimha  model  simply  expresses 
the  intermittency  factor.  If,  which  indicates  the  probability  of  the  flow 
being  turbulent  at  some  point  as 


If(x)  =  1  -  exp 


(28) 
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where  0  is  a  constant  equal  to  0. 412  and  X  is  a  measure  of  the  extent 
of  the  transition  region  characterized  by 

x  =  *]If  =  0.75  "  X]If  =  0.25  (29) 

The  beginning  of  the  transition  zone,  where  lf(0)  =  0,  was  estab¬ 
lished  by  Dhawan  and  Narasimha  (Ref.  18)  as  being  located  at  the 
transition  onset  distance  x  =  Xj.  In  order  to  describe  the  extent  of  the 
transition  zone,  Masaki  and  Yakura  (Ref.  1)  assume  that  the  end  of 
transition  is  achieved  at  x  =  X»p  where 

If(XT)  *  0.97  (30) 

Substituting  this  assumption  into  Eq.  (28)  gives  the  length  of  the  tran¬ 
sition  zone  as 


AX  =  XT  -  Xt  =  2.96  x 


(31) 


In  order  to  evaluate  x,  recall  that  Eq.  (27)  mathematically  defines  the 
end  of  transition  in  terms  of  the  transition  onset  distance  Xj-  as 

Ree,Xt  “  2Ree,Xt  (32) 


which  implies  physically  that  the  transition  zone  length  Ax  is  equal  to 
the  transition  onset  distance  X^  for  the  special  case  of  supersonic  or 
hypersonic  flow  over  a  sharp  cone  or  sharp  flat  plate  where  the  inviscid 
local  flow  parameters  (p,  Te,  Pe,  Ue)  remain  constant  over  the  entire 
body.  Under  this  restriction 

Ax  =  Xt  (33) 

so  that  X  becomes,  from  Eq.  (31) 

Xt 

x  "  “OF  <34> 


With  the  above  restrictions  the  intermittency  factor  given  by  Eq.  (28) 
takes  the  final  form 


If(x)  =  1  -  exp 


-  0.412  (2.96)2 


(35) 


where 


Xt 
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Although  the  original  derivation  of  the  intermittency  factor  was 
based  on  experiments  performed  in  low  speed  incompressible  flow, 

Dhawan  and  Narasimha  {Ref.  18)  establish  its  validity  for  use  in  pre¬ 
diction  of  supersonic  skin  friction  on  a  flat  plate.  To  the  author's 
knowledge,  no  work  has  been  done  on  the  effects  of  pressure  gradients 
on  the  If(x)  distributions  in  either  incompressible  or  compressible  flows, 
although  Dhawan  and  Narasimha  speculate  that  pressure  gradients  would 
chiefly  influence  the  manner  of  spot  growth  in  the  initial  period  near  If  =  0. 
Another  point  of  interest  is  the  distribution  of  If  across  the  boundary 
layer,  viz,  normal  to  the  surface.  Dhawan  and  Narasimha  state  that 
although  this  If (y)  variation  is  probably  of  importance  to  the  detailed 
structure  of  the  turbulent  motion,  the  value  of  If  near  the  wall  is  the 
characteristic  property  of  importance  for  the  transition  region  and  the 
If(y)  variation  has  only  a  secondary  influence  in  determining  the  mean 
velocity  profiles  in  transition.  Hence,  for  the  hypersonic  sharp  cone 
flows  of  interest  in  the  present  investigation,  Eq.  (35)  is  used  for  the 
intermittency  factor  distribution  in  order  to  assess  its  validity  for  use 
in  such  regimes. 

In  the  present  work  the  intermittency  factor  multiplies  the  eddy  vis¬ 
cosity  in  the  governing  boundary -layer  equations  (5)  and  (6).  This  appli¬ 
cation  follows  Batchelor  (Ref.  19)  who  found  that  this  type  of  analytical 
treatment  yielded  the  best  agreement  with  experiments  for  the  case  of 
two-dimensional  incompressible  far  wake  flows  where  one  is  concerned 
with  variations  across  the  wake.  The  physical  justification  for  using  the 
intermittency  factor  in  this  manner  lies  in  the  basic  physical  definition 
of  the  intermittency  factor  as  representing  the  fraction  of  time  any  point 
spends  in  turbulent  flow,  i.  e. ,  the  probability  of  the  flow  being  turbulent 
at  a  given  point.  If  If  =  0,  the  flow  is  completely  laminar,  but  if  If  =  1, 
the  flow  is  fully  turbulent.  The  present  treatment  of  the  intermittency 
factor-eddy  viscosity  product  in  the  governing  equations  simply  means 
that  the  (fully  turbulent)  eddy  viscosity  is  multiplied  by  a  damping  factor 
(the  intermittency  factor)  in  order  to  achieve  a  smooth  and  continuous 
variation  from  laminar  to  turbulent  flow  through  the  transition  region.  A 
similar  approach  has  been  used  by  Martellucci,  Rie,  and  Sontowski 
(Ref.  6)  who  varied  the  eddy  viscosity  linearily  with  respect  to  distance 
through  the  transition  zone.  In  a  sense  this  may  be  considered  an  inter¬ 
mittency  factor  type  of  treatment  where  the  intermittency  factor  is  a 
linear  function  of  distance.  However,  Martellucci,  Rie,  and  Sontowski 
did  not  attempt  to  assess  the  validity  of  this  approach;  no  results  were 
presented  of  transition  zone  predictions  relative  to  experimental  data. 

In  this  respect  the  present  work  is  the  first,  to  the  author's  knowledge, 
concerning  analytical  techniques  applicable  to  the  transition  region  where 
the  governing  boundary-layer  equations  are  numerically  integrated  on  a 
digital  computer  yielding  predictions  which  are  then  compared  with  ex¬ 
perimental  data  in  order  to  establish  their  validity. 
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2.5  COORDINATE  TRANSFORMATION 

In  order  to  facilitate  numerical  integration  of  the  governing  boundary 
layer  equations  (4),  (5),  and  (6),  it  is  convenient  to  transform  them  to  a 
coordinate  system  that  removes  the  singularity  at  x  =  0  and  stretches  the 
coordinate  normal  to  the  flow  direction,  as  is  usually  done  in  laminar 
flow.  The  coordinate  transformation  used  in  the  present  work  is  the  well 
known  Lees -Dorodnitsyn  (Ref.  9,  p.  31)  transformation.  The  new  inde¬ 
pendent  variables  introduced  are 

X 

5(x)  =  f  peye  b'e  rb2j  dx 

Jo 

(36) 

,  .  PeW  f  £_  . 
n(x’y)  a  m  ^ 

(37) 

and  hence  derivatives  become 

!_=  p  v  U  n2^  i-  +  ML 

3x  peye  e  T)  35  ax  3n 

(38) 

a  p  3 

ay  "  JzT  3n 

(39) 

Introducing  a  stream  function  ip(x,  y)  in  such  a  manner  as 
satisfy  the  continuity  equation  (4) 

to  identically 

»  “V1=!y 

(40) 

P  V  rb°='3'x 

(41) 

and  defining  a  transformed  stream  function  f(|,  17)  such  that 

¥  (x,y)=  /2T  fU,n) 

(42) 

the  governing  boundary -layer  equations  (Eqs.  (5)  and  (6))  become,  in 
the  transformed  f,rj  plane: 
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MOMENTUM 


U*  f"  )  v  ff "  +  B[e  -  (f')2] 


=  2c 


r  .  f~  lf- 
T  ac  T  ac 


(43) 


ENERGY 


a**  s'  '  +  fe'  +  (y-i)mc2  [i*  (f"  )2 


-  eef'J  =  2c 


ii  _  11 
1  as  as 


e'  + 


with  the  new  dependent  variables 


ef1  dTo 

T«  3T 


(44) 


f'  (S,n)  =  |j- 


e(s,n)  «  r- 
e 


(45) 

(46) 


where  primes  denote  differentiation  with  respect  to  the  independent 
variable  r],  i.  e. , 

_  a2f 

17 


0'  =  ^ 
6  an 


The  following  definitions  apply  to  the  above  equations 


l*  = 


PU 

peve 


1  +  If  S’ 


(47) 


„**  =  PiL_  |i  +  i  L  PL_| 

Peue  l  f  p  Prt  I 


(48) 
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2?  dUe 

B  “  U£  '  (49) 

_Ue _ 

e  ^7  (50) 

where  the  inviscid  momentum  equation 


has  been  used  to  relate  the  inviscid  velocity  Ue  to  the  imposed  inviscid 
static  pressure. 

The  physical  boundary  conditions  given  by  Eqs.  (9)  through  (13) 
become,  in  the  transformed  variables: 


MOMENTUM 


fU.Ti  =  0)  =  o 

(52) 

f *(€,n  -  0)  -  0 

(53) 

lim  f'(?,n)  =  1 
n  -*■  00 

(54) 

ENERGY 

0(5, n  =  0)  =  ew  U) 

(55) 

lim  0(5, n)  =  1 

p  -V  00 

(56) 

The  local  convective  heat  flux  at  the  body  surface  (y 
by  the  well-known  Fourier  law 

=  0)  is  given 

11  = '  ^  8ylw 

(57) 

which  becomes,  in  the  transformed  coordinates. 

q  =  1, CpTe  0,{5jTj  =  Q) 

Prw  JTi 

(58) 
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where 


“w 


Pw^w 

peMe 


(59) 


In  a  similar  manner,  the  local  wall  shearing  stress  may  be  written  in 
terms  of  the  transformed  coordinates  as 


lw 


% 


3  U 

3Ww 


(60) 

=  =  0) 

Various  nondimensional  quantities  of  interest  are  defined  as  follows: 


LOCAL  REYNOLDS  NUMBER  BASED  ON  FREE-STREAM  CONDITIONS 

„  PcUooX 

Re  Y  *  - 

oo, X  p 


(61) 


LOCAL  REYNOLDS  NUMBER  BASED  ON  INVISCID  EDGE  CONDITIONS 

P@U@  x 


Re 


e,x 


(62) 


LOCAL  SKIN  FRICTION  COEFFICIENT  REFERENCED  TO  FREE-STREAM 
CONDITIONS 

TW 


Cf  =  T  2 
“  T  P_  U  2 


(63) 


LOCAL  STANTON  NUMBER  BASED  ON  INVISCID  EDGE  CONDITIONS  AND 
ADIABATIC  WALL  TEMPERATURE 


St, 


e ,  aw  Peue^p  ^aw  ~  ^ 


(64) 
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where 


rf  +  j2 —  (1  -  rf) 
'  '0,«  T 


(65) 


with  rf  the  recovery  factor  which  must  be  defined  relative  to  the  flow 
under  examination. 


2.6  NUMERICAL  SOLUTION 

In  the  present  work,  numerical  solution  of  the  simultaneous,  non¬ 
linear,  parabolic,  partial  differential  equations  (Eqs.  (43)  and  (44))  is 
performed  by  obtaining  linear  finite -difference  equivalents  of  the  equa¬ 
tions  and  solving  these  using  an  iterative,  marching,  tridiagonal 
matrix  method.  Details  of  the  finite -difference  procedure  are  given  in 
Appendix  III. 

The  digital  computer  code  is  written  in  FORTRAN  63  for  use  on  a 
CDC  1604-B  machine.  Solution  time  including  printout  is  approximately 
20  minutes  on  the  1604-B  machine  for  a  body  divided  into  100  stations. 
No  numerical  stability  problems  have  been  encountered  with  the  present 
finite-difference  approach  because  of  its  implicit  nature. 

Much  debate  currently  exists  about  the  relative  merits  of  different 
numerical  techniques  for  solution  of  the  turbulent  boundary-layer  equa¬ 
tions;  a  good  example  may  be  found  in  Paper  21  of  Ref.  20  which  is  a 
written  transcription  of  the  verbal  panel  and  general  discussion  portion 
of  a  1968  symposium  on  compressible  turbulent  boundary  layers. 
General  reviews  of  this  subject  have  been  given  by  Spalding  (Ref.  21), 
McDonald  (Ref.  22),  Beckwith  (Ref.  23),  and  Smith  (Ref.  24). 


SECTION  III 

COMPARISON  OF  PRESENT  NUMERICAL  RESULTS  WITH  EXPERIMENTAL  DATA 


As  pointed  out  by  Smith  (Ref.  24)  in  his  recent  review  of  boundary- 
layer  research  over  the  past  decade,  the  ultimate  test  of  any  theoretical 
analysis  of  the  turbulent  boundary  layer  is  the  end  accuracy  (and  not  the 
intrinsic  logic)  since  empiricism  is  involved.  Hence,  with  this  in  mind, 
comparison  of  experimental  data  with  the  present  eddy  viscosity- 
intermittency  factor  approach  to  prediction  of  transitional  heating  under 
hypersonic  conditions  is  mandatory.  Since  the  current  effort  is  limited 
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solely  to  sharp  cone  geometry,  the  following  additional  restrictions  are 
imposed  at  this  point  in  order  to  clearly  define  the  experimental  flow 
field  for  comparison  with  the  analytical  model: 

1.  Hypersonic  flow  at  zero  angle  of  attack. 

2.  Uniform  free  stream,  i.  e.  ,  no  source  flow  effects. 

3.  Cold  wall  at  constant  uniform  temperature. 

4.  Transition  zone  well  defined  by  heat -transfer  data. 

Under  the  above  restrictions,  a  review  of  unclassified,  published 
literature  on  hypersonic  sharp  cone  experimental  investigations  re¬ 
vealed  the  six  acceptable  sources  (Refs.  25  through  30)  listed  in  Table  I, 
Appendix  II.  Although  this  list  appears  to  be  very  meager,  the  majority 
of  published  investigations  do  not  meet  the  rather  stringent  restrictions 
imposed  above  and  hence  must  be  eliminated  on  this  basis.  For  example, 
the  work  by  Nagamatsu  and  Sheer  (Ref.  31)  contains  source  flow  effects, 
and  the  wall  temperature  is  not  uniform  in  the  study  by  Everhart  and 
Hamilton  (Ref.  32);  Stetson  and  Rushton  (Ref.  33)  present  no  usable  heat- 
transfer  data  defining  transition.  By  only  using  experimental  data  from 
Table  I,  meaningful  comparisons  between  theory  and  experiment  can  be 
made  with  some  confidence  since  the  assumptions  under  which  the  ana¬ 
lytical  model  is  formulated  match  the  restrictions  placed  on  acceptable 
experimental  data. 

Table  II  lists  the  inviscid  outer-edge  conditions  for  input  to  the 
boundary -layer  analysis  corresponding  to  the  various  flows  in  Table  I. 
The  inviscid  conditions  are  determined  from  numerical  solution  of  the 
inviscid  conical  flow  equations  following  Sims  (Ref.  34)  for  a  given 
sharp  cone  half-angle  and  free-stream  Mach  number.  Consistent  with 
the  gas  model  adopted  in  Section  2.1,  a  perfect  gas  with  y  =  1. 40  is 
used  in  the  inviscid  analysis. 

Of  the  data  in  Table  I,  the  investigation  by  McCauley,  Saydah,  and 
Bueche  (Ref.  25)  is  one  of  the  best  documented  with  respect  to  test  condi¬ 
tions  and  analysis  of  data.  The  test  was  conducted  in  the  Hypersonic 
Wind  Tunnel  (C)  of  the  von  Karman  Gas  Dynamics  Facility,  Arnold  Engi¬ 
neering  Development  Center,  at  a  nominal  Mach  number  of  10.  Test 
conditions  are  given  in  Table  I.  The  objective  of  the  experimental 
investigation  was  to  determine  the  effects  of  various  diameter  spherical 
trips  on  boundary -layer  transition  along  sharp  and  spherically  blunted 
cones.  Figure  2  shows  a  comparison  of  results  from  the  present  theory 
relative  to  experimental  heat -transfer  measurements  for  both  natural 
transition  and  roughness -induced  transition.  In  general,  the  agreement 
between  experimental  data  and  results  from  the  present  theory  is  good. 
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The  transition  region  is  well  defined  by  the  present  theory  which,  it 
should  be  recalled,  requires  that  the  location  of  onset  to  transition  be 
prescribed  as  shown  in  Fig,  2.  Both  sets  of  experimental  data  pre¬ 
sented  in  Fig.  2  clearly  show  a  factor  of  about  2  in  the  ratio  X-p/Xj.  as 
previously  discussed  in  Section  2.4  and  used  in  the  present  analytical 
model.  For  the  tripped  condition,  note  that  fully  turbulent  flow  exists 
over  approximately  half  the  cone  but,  for  the  natural  transition  condi¬ 
tion,  the  flow  remains  transitional  over  the  entire  aft  half  of  the  body. 

It  is  the  large  extent  of  transitional  flow  in  the  latter  case  which  poses 
the  thermal  design  problem  for  lifting  entry  vehicles. 

Since  the  agreement  between  results  from  the  present  theory  and 
experiment  as  shown  in  Fig.  2  indicates  that  the  present  theory  is  indeed 
applicable  to  prediction  of  transitional  heating  on  a  sharp  cone  in  hyper¬ 
sonic  flow,  it  is  now  in  order  to  examine  certain  details  of  the  flow  field. 
Figures  3  through  7  present  calculated  velocity,  static  temperature, 

Mach  number,  stagnation  temperature,  and  pitot  pressure  profiles  at 
three  locations  along  the  sharp  cone  for  the  roughness -induced  transition 
condition  where  the  onset  of  transition  occurs  at  the  body  location 
x/L  =  0.  24.  This  choice  of  flow  conditions  has  been  made  to  permit 
presentation  of  laminar  (x/L  =  0.  20),  transitional  (x/L  =  0.40),  and 
fully  turbulent  (x/L  =  0.  80)  profiles;  the  station  x/L  =  0.  40  for  the  tran¬ 
sitional  profile  corresponds  to  approximately  the  location  of  calculated 
maximum  heating  in  Fig.  2.  The  different  character  of  the  laminar  and 
turbulent  profiles  is  apparent.  The  transitional  profiles  exhibit  a  larger 
gradient  in  flow  properties  near  the  wall  than  do  the  fully  turbulent  pro¬ 
files.  This  phenomenon  is  discussed  in  more  detail  later. 

In  order  to  present  the  differences  between  laminar  profiles  and 
fully  turbulent  profiles  at  a  given  body  location.  Figs.  8  and  9  show 
calculated  velocity  and  static  temperature  distributions  for  the  two  types 
of  flow  at  x/L  =  0.  80.  Figure  8  reveals  the  much  fuller  turbulent  profile 
relative  to  the  laminar  profile  as  expected,  whereas  Fig.  9  shows  that  the 
peak  static  temperature  in  the  laminar  boundary  layer  occurs  much. far¬ 
ther  from  the  wall  than  the  peak  turbulent  static  temperature,  about  a 
factor  of  3  in  distance  normal  to  the  surface.  This  factor  of  3  in  loca¬ 
tion  of  peak  static  temperature  is  somewhat  reflected  in  the  thickness 
of  the  turbulent  boundary  layer  being  approximately  twice  the  thickness 
of  the  laminar  boundary  layer  as  can  be  seen  by  close  inspection  of 
Figs.  8  and  9. 

Another  comparison  of  interest  relative  to  laminar  and  transitional 
flow  profiles  near  the  onset  of  transition  is  presented  in  Figs.  10  and  11 
for  the  station  x/L  =  0.  30.  The  .transition  process  to  turbulent  flow 
affects  the  entire  profile.  This  observation  is  in  quantitative  agreement 
with  the  experimental  observations  of  Maddalon  and  Henderson  (Ref.  35) 
as  well  as  Potter  and  Whitfield  (Ref.  36)  for  natural  transition. 
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In  order  to  understand  the  profile  behavior  discussed  above,  one 
must  keep  in  mind  that  the  eddy  viscosity-intermittency  factor  product 
governs  the  transitional  and  turbulent  flow  calculations  in  the  current 
approach.  Figure  12  gives  the  intermittency  factor  variation  over  the 
body  as  governed  by  Eq.  (35),  and  Fig.  13  gives  the  calculated  eddy 
viscosity-intermittency  factor  product  distribution  across  the  boundary 
layer  at  selected  stations.  As  is  apparent  from  Fig.  13,  the  eddy  vis¬ 
cosity  reaches  its  maximum  value  in  the  outer  region  of  the  boundary 
layer  with  e  »  n  for  fully  turbulent  flow.  Note  that  for  sufficiently 
large  values  of  x/L  such  that  If  »  1,  the  ratio  €  In  at  a  given  y/L  value 
in  the  wall  region  approaches  a  common  value  independent  of  the  x  loca¬ 
tion  as  may  be  seen  by  comparing  the  curves  for  x/L  =  0.  40  and  0. 80 
in  Fig.  13  at  values  of  y/L  <  5.0  x  10"4.  In  the  transition  region  where 
If  <  1,  the  value  of  the  e  In  ratio  near  the  wall  is  approximately  the  fully 
turbulent  value  at  the  same  y/L  location  scaled  down  by  the  If  ratio. 

For  example,  at  y/L  =  5.  0  x  10“^,  e  In  =0.  14  for  x/L  =  0.  30  and 
e  In  =0.  68  for  x/L  =  0.  80,  a'factor  of  approximately  five  difference. 
From  Fig.  12,  If  =  0.  20  at  x/L  =  0.  30  and  If  =  1.0  at  x/L  =  0.  80  which 
reveals  the  factor  of  five  scaling  noted  above.  Thus  the  effect  of  the 
intermittency  factor  in  the  wall  region  is  simply  a  scaling  of  the  e  In 
ratio.  In  the  outer  regions  of  the  boundary  layer  where  the  change  has 
been  made  from  the  inner  to  the  outer  expression  for  calculating  the 
eddy  viscosity,  scaling  of  the  c  In  ratio  is  not  applicable  as  can  be  seen 
by  reference  to  Fig.  13.  Hence  the  character  of  the  transitional  zone 
description  in  the  present  analysis  as  presented  in  Figs.  10  and  11  can 
be  directly  traced  to  the  treatment  of  the  eddy  viscosity-intermittency 
factor  product  discussed  above. 

One  well -documented  phenomenon  in  transitional  bound  ary -layer 
flows  is  the  experimentally  observed  overshoot  in  peak  heat-transfer 
rate  near  the  end  of  transition  relative  to  fully  turbulent  boundary- 
layer  theory.  Figure  14  shows  a  comparison  of  results  from  the  pres¬ 
ent  method  of  analysis  using  the  eddy  viscosity-intermittency  factor 
treatment  of  transition  and  results  using  an  instantaneous  step  transi¬ 
tion  from  laminar  to  turbulent  flow  at  a  given  location.  The  present 
method  clearly  shows  an  overshoot  in  heat-transfer  rate.  This  over¬ 
shoot  in  heat -transfer  rate  is  attributable  to  the  transitional  profile 
behavior  discussed  previously  in  connection  with  Figs.  3  through  7. 
Similar  comments  apply  to  the  wall  skin  friction  coefficient  as  shown  by 
Fig.  15.  One  important  point  to  note  from  Fig.  14  is  that  the  calculated 
heat-transfer  rates  relax  to  common  downstream  values  irrespective  of 
the  transition  zone  treatment;  the  same  is  true  for  the  wall  skin  friction 
coefficients  as  can  be  seen  from  Fig.  15.  This  behavior  is  to  be 
expected  because  of  the  parabolic  mathematical  character  of  the  govern¬ 
ing  boundary -layer  equations.  Further  note,  from  Fig.  16,  that  the 
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instantaneous  step  transition  treatment  gives  a  larger  value  for  the 
boundary -layer  thickness  downstream  than  the  present  intermittency 
factor-eddy  viscosity  approach.  Sufficiently  far  downstream  both 
models  give  the  same  value. 

One  method  often  used  by  experimentalists  to  detect  boundary-layer 
transition  is  a  pitot  tube  mounted  close  to  the  surface  of  the  body.  Fig¬ 
ure  17  shows  the  calculated  pitot  pressure  variation  along  the  body  at 
two  fixed  locations  normal  to  the  surface.  The  variation  through  the 
transition  zone  is  well  defined,  and  the  general  trends  are  similar  to 
those  for  the  distributions  of  heat-transfer  rate  and  skin  friction  coeffi¬ 
cient  from  Figs.  14  and  15.  Although  the  distance  y/L  =  4.  0  x  10"^ 
may  not  be  physically  realistic  from  the  standpoint  of  actually  construct¬ 
ing  a  working  probe,  the  above  results  indicate  that  pitot  pressure  meas¬ 
urements  do  indeed  offer  a  valid  means  of  detecting  transition  under 
hypersonic  cold  wall  conditions. 


Another  technique  often  used  by  the  experimentalist  is  the  measure¬ 
ment  of  pitot  pressure  profiles  across  the  boundary  layer,  which  are 
used  in  conjunction  with  the  Crocco  relation  (Ref.  9,  pp.  33  and  182) 


(66) 


to  infer  velocity  profiles.  The  Crocco  relation  (Eq.  (66))  represents  a 
solution  to  the  governing  laminar  and  turbulent  boundary -layer  equations 
for  the  special  case  of  zero  pressure  gradient,  constant  wall  temperature, 
and  both  laminar  and  turbulent  Prandtl  numbers  equal  to  unity.  These 
conditions  are  satisfied  in  the  present  analysis  of  hypersonic  sharp  cone 
flows  (except  for  the  Prandtl  number  restriction).  Figure  18  presents  a 
comparison  of  the  present  theory  with  the  Crocco  relation  written  in 
terms  of  the  present  transformed  variables  as 


9  =  9W  +  f'  0  -  9W) 


(67) 


where 


9 
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(68) 

(69) 


with  f'  =  u/Ue  from  Eq.  (45).  In  general  the  agreement  is  excellent 
which  indicates  that  the  Prandtl  number  and  the  turbulent  Prandtl  num 
ber  have  only  a  small  influence  on  the  solution  when  expressed  in  this 
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form.  It  should  be  pointed  out,  however,  that  this  observation  is  being 
made  for  cold  wall  conditions  only. 

The  investigations  by  Wallace  (Refs.  37  and  38)  are  an  excellent 
example  of  the  approach  given  in  the  previous  paragraph.  In  order  to 
assess  whether  calculated  velocity  profiles  were  laminar  or  turbulent, 
Wallace  compared  them  with  power  law  profiles  of  the  form 

(d)  (70) 

since  turbulent  boundary  -layer  profiles  are  usually  associated  with 
small  values  of  n,  typically  n  =  1/7  to  1/10.  Figure  19  presents  repre¬ 
sentative  power-law  profiles  relative  to  the  calculated  velocity  profile 
from  the  present  theory  at  a  body  location  x/L  =0.  80  where  a  fully 
turbulent  boundary  layer  is  present.  Above  y/6  =  0.20,  the  present 
theory  is  seen  to  agree  closely  with  the  n  =  1/5  power  law.  It  is 
interesting  to  note  in  this  respect  that  Pinckney  (Ref.  39)  found  that 
a  1/5  power  law  velocity  profile  yielded  good  agreement  with  experi¬ 
mental  static  temperature -velocity  profiles  for  flat  plate  compressible 
turbulent  boundary  layers  with  heat  transfer  in  the  free -stream  Mach 
number  range  from  0.  85  to  8.  18. 

So  far  the  present  work  has  been  devoted  to  an  analysis  of  one 
particular  sharp  cone  flow  under  hypersonic  conditions.  Returning  to 
the  experimental  data  in  Table  I,  it  is  now  in  order  to  compare  results 
from  the  present  theory  with  the  remaining  data.  In  all  cases  the  com¬ 
parisons  will  be  made  in  terms  of  the  same  parameters  given  in  the 
references  so  as  to  avoid  error  in  reduction  of  data  to  a  different  form. 

Shown  in  Fig.  20  is  a  comparison  of  a  heat -transfer  rate  parameter 
calculated  using  the  present  theory  with  the  data  presented  by  DiCristina 
(Ref.  26)  for  Mach  10.  2  hypersonic  flow.  This  test  was  conducted  in  the 
Hypersonic  Wind  Tunnel  (C)  of  the  von  Karman  Gas  Dynamics  Facility, 
Arnold  Engineering  Development  Center;  test  conditions  are  given  in 
Table  I  with  corresponding  inviscid  edge  conditions  in  Table  H.  Agree¬ 
ment  of  results  from  the  present  theory  and  experiment  is  excellent.  A 
transition  zone  length  equal  to  Xj.  is  clearly  seen  in  these  data  obtained 
under  conditions  for  which  natural  transition  occurred. 

Figure  21  presents  the  experimental  data  of  Sanator,  DeCarlo,  and 
Torrillo  (Ref.  27)  and  results  from  the  present  theory.  This  test  was 
performed  in  the  3 6 -in.  hypersonic  wind  tunnel  at  the  Republic  Aviation 
Corporation;  test  conditions  are  given  in  Table  I  with  corresponding 
inviscid  edge  conditions  in  Table  II.  Spherical  roughness  elements  were 
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used  to  trip  the  boundary  layer  in  this  test,  with  the  location  of  these 
elements  shown  in  Fig.  21.  Again,  the  agreement  of  results  from  the 
present  theory  and  experiment  is  excellent  with  fully  turbulent  flow 
being  established  over  approximately  half  the  cone.  The  transitional 
zone  appears  well  defined  by  the  present  theory. 

A  very  extensive  program  for  measurement  of  hypersonic  turbulent 
heat  transfer  on  a  highly  cooled  sharp  cone  has  been  under  way  for  several 
years  at  the  United  States  Naval  Ordnance  Laboratory  using  their  Mach  5 
hypersonic  wind  tunnel.  The  results  of  these  investigations  have  N»en 
published  in  reports  by  Wilson  and  Fisher  (Ref.  28)  and  Wilson  (R 
Figure  22  shows  a  comparison  of  results  from  the  present  theory  with 
experimental  data  from  Fig.  6a  of  Ref.  28  and  Fig.  13c  of  Ref.  29;  test 
conditions  and  inviscid  outer  edge  conditions  are  given  in  Tables  I 
and  II,  respectively.  As  shown  in  Table  I,  the  test  conditions  are  almost 
identical  for  the  two  sets  of  data.  The  results  from  the  present  theory 
agree  best  with  the  data  from  Ref.  2  9  in  the  fully  turbulent  region  whereas 
the  data  from  Ref.  28  appear  to  best  define  the  beginning  of  the  transition 
region.  Since  the  two  sets  of  data  are  not  mutually  consistent  it  is  diffi¬ 
cult  to  assess  if  indeed  the  present  theory  is  overestimating  the  fully 
turbulent  heat-transfer  rate  for  these  flow  conditions.  Comparisons  not 
presented  herein  of  results  from  the  present  theory  with  other  experi¬ 
mental  data  in  Refs.  28  and  29  for  different  flow  conditions  reveal  that  in 
all  cases  the  present  theory  gives  values  higher  than  the  experimental 
data.  It  should  be  noted  that  Mayne  and  Dyer  (Ref.  7)  also  obtained 
higher  values  than  the  fully  turbulent  experimental  data  from  Ref.  29  as 
can  be  seen  by  reference  to  Fig.  6  of  their  paper. 

Another  experimental  investigation  at  a  low  hypersonic  Mach  number 
(4.  95)  is  reported  by  Julius-(Ref.  30).  This  test  was  conducted  in  the 
9-in.-diam  blowdown  axisymmetric  jet  at  the  NASA  Langley  Research 
Center;  test  conditions  are  given  in  Table  I.  The  inviscid  edge  conditions 
given  in  Table  II  for  this  flow  condition  are  based  upon  an  average  of  the 
experimental  pressure  distribution  along  the  body  shown  in  Fig.  3  of 
Ref.  30.  Because  of  tunnel  effects  the  model  experienced  a  small  de¬ 
crease  in  surface  pressure  over  the  aft  third  of  the  body.  The  compari¬ 
son  shown  in  Fig.  23  reveals  excellent  agreement  of  results  from  the 
present  theory  and  experiment  up  to  the  point  where  the  (favorable)  pres¬ 
sure  gradient  tends  to  reduce  the  experimental  turbulent  heat -transfer 
rate  to  a  value  below  the  theoretical  prediction.  In  the  laminar  region 
the  experimental  data  appear  high  relative  to  results  from  the  present 
theory.  This  is  caused  by  the  experimental  pressure  in  the  laminar 
region  being  higher  than  the  averaged  pressure  used  in  the  present  calcu¬ 
lations. 
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SECTION  IV 

CONCLUDING  SUMMARY 


In  the  above  comparisons  of  experimental  results  with  numerical 
calculations  from  the  present  theory,  excellent  agreement  was  obtained 
in  most  cases.  The  present  approach  appears  to  yield  a  reasonably 
accurate  description  of  transitional  heating  on  slender  sharp  cones  in 
the  Mach  number  range  from  5  to  10  under  cold  wall  conditions.  It  is 
important  to  keep  in  mind  that  application  of  the  present  theory  requires 
a  priori  specification  of  the  onset  to  transition  location,  X^,  along  the  body. 
In  the  work  herein,  this  location  was  determined  based  upon  examination 
of  the  particular  experimental  data  under  investigation.  The  factor  of  2 
used  to  define  the  end  of  transition  relative  to  the  transition  onset  dis¬ 
tance  appears  to  be  adequate  for  description  of  sharp  cone  transitional 
heat  transfer  under  cold  wall  conditions  in  the  hypersonic  flow  regime. 

These  favorable  results  indicate  that  the  intermittency  factor -eddy 
viscosity  approach  coupled  with  accurate  numerical  solution  of  the 
governing  turbulent  boundary-layer  equations  may  become  a  powerful 
tool  in  analysis  techniques  for  transitional  heating  on  lifting  body  con¬ 
figurations.  Much  work  remains  to  be  done  in  this  area  as  will  be  dis¬ 
cussed  in  the  next  section. 


SECTION  V 
FUTURE  STUDIES 


The  continuing  need  for  understanding  turbulent  boundary  layers  is 
more  experiments  of  a  so-called  microscopic  nature,  i.  e. ,  detailed 
measurements  of  the  flow  field  structure  performed  in  parallel  with 
accurate  application  of  the  theories  now  in  existence.  Above  Mach  num¬ 
ber  5  acceptable  detailed  measurements  of  the  turbulent  boundary-layer 
structure  are  almost  nonexistent  for  flows  over  bodies  of  aerodynamic 
interest.  Much  evidence  (see  Ref.  40)  seems  to  indicate,  at  least  in 
the  high  supersonic  or  low  hypersonic  Mach  number  range,  that  the  gen¬ 
eral  concept  of  a  turbulent  mechanism  (such  as  used  in  the  present  work) 
unaffected  by  the  supersonic  nature  of  the  mean  stream  is  acceptable. 
The  characteristics  of  turbulence  in  both  incompressible  and  compress¬ 
ible  flow  seem  to  depend  only  on  the  kinematics  of  the  flow.  Hence  it 
may  well  be  that  one  can  safely  apply  the  present  simple  structural 
hypothesis  with  confidence  in  the  hypersonic  regime;  at  least  the  present 
work  indicates  that  an  eddy  viscosity  approach  yields  accurate  heat- 
transfer  predictions  for  hypersonic  sharp  cone  transitional  and  turbulent 
boundary -layer  flows. 
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With  respect  to  extension  of  the  present  eddy  viscosity- intermittency 
factor  treatment  of  transition,  it  is  the  author's  opinion  that  the  following 
investigations  would  be  of  great  value  in  establishing  its  validity  under 
flow  conditions  different  from  the  current  investigation: 

1.  Comparison  of  theory  and  experiment  with  respect  to 
heat  transfer  and  skin  friction  on  flat  plates  in  super-, 
sonic  and  hypersonic  flow. 

2.  Comparison  of  theory  and  experiment  with  respect  to 
heat  transfer  and  skin  friction  on  blunt  nosed  bodies 
in  supersonic  and  hypersonic  flow. 

3.  Comparison  of  theory  and  experiment  with  respect  to 
the  recovery  factor  under  adiabatic  wall  conditions  in 
both  supersonic  and  hypersonic  flow. 

Again  detailed  measurements  of  the  turbulent  boundary -layer  structure 
in  each  of  the  above  suggested  investigations  are  direly  needed  in  order 
to  strictly  assess  the  theoretical  model.  Another  important  result  of 
the  above  suggested  investigations  would  be  the  compilation  of  sufficient 
experimental  data  for  the  body  and  flow  under  consideration  to  permit 
determination  of  accurate  transition  zone  descriptive  models.  There  is 
no  reason  to  believe  that  the  factor  of  2  transition  zone  model  used  in 
the  present  sharp  cone  investigation  will  remain  valid  for,  say,  a  curved 
flat  plate  having  an  adiabatic  wall  under  low  supersonic  conditions.  It 
may  well  be  that  additional  experiments  must  be  designed  in  order  to 
carefully  study  the  transition  zone  under  controlled;flow  conditions  for 
various  body  geometries.  In  any  respect,  comparison  of  theory  with 
experiments  is  the  key  to  understanding  the  applicability  of  the  current 
eddy  viscosity-intermittency  factor  prediction  technique. 

Another  area  worthy  of  future  work  is  the  application  of  turbulent 
transport  theory  to  numerical  calculation  of  transitional  boundary -layer 
structure.  A  recent' paper  by  Bradshaw,  Ferriss,  and  Atwell  (Ref.  41) 
successfully  applies  the  turbulent  energy  equation  (which  is  an  equation 
for  the  rate  of  change  of  turbulent  intensity  along  a  streamline)  to 
incompressible  turbulent  flows  both  with  and  without  pressure  gradients. 
While  this  method  employs  several  empirical  relations,  it  is  more  gen¬ 
eral  than  the  eddy  viscosity  boundary -layer  approach.  Prediction  of 
transition  should  conceptually  be  possible  using  this  technique.  The 
only  attempt  published  to  date  to  solve  the  compressible  boundary  -layer 
problem  by  this  approach  is  that  of  Bradshaw  and  Ferriss  (Ref.  42)  who 
extended  their  method  for  incompressible  flow  directly  to  adiabatic  flows 
with  M0  <  4.  It  certainly  would  be  of  interest  to  ascertain  if  the  present 
intermittency  factor -eddy  viscosity  approach  is  contained  within  the  basic 
framework  of  turbulent  transport  theory. 
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The  present  analysis  reported  herein  covers  the  first  phase  of  a 
general  analytical  investigation  having  as  the  end  objective  accurate 
numerical  prediction  techniques  for  laminar,  transitional,  and  turbu¬ 
lent  boundary  layers  in  hypersonic  flow  along  the  most  windward  stream¬ 
line  of  arbitrarily  shaped  bodies  at  angle  of  attack.  Subsequent  reports 
will  document  the  findings  of  this  ongoing  work. 
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Fig.  2  Comparison  of  Present  Results  with  Experimental  Data 
of  McCauley,  Saydah,  and  Bueche  (Ref.  25) 
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Fig.  3  Calculated  Velocity  Profiles 
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Fig.  8  Comparison  of  Laminar  and  Fully  Turbulent  Velocity  Profiles 
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Fig.  11  Comparison  of  Laminar  and  Transitional  Static  Temperature  Profiles 
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Fig.  13  Eddy  Viscosity  Distributions  at  Various  Body  Locations 
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Fig.  14  Effects  of  Various  Transition  Zone  Models  on  Heat-T ransfer- Rate  Distribution 


48 


AEDC-TR-70-210 


6.0-deg  Half-Angle  Sharp  Con 
M  -10.17,  Re  ffi-1.6xltf 

VW0'29'  L-47'38in- 


Fig.  15  Effects  of  Various  Transitional  Zone  Models  on  Skin  Friction 
Coefficient  Distribution 
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Fig.  16  Effects  of  Various  Transition  Zona  Models  on  Boundary-Layer 
Thickness 
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Fig.  21  Comparison  of  Present  Results  with  Experimental  Data  of  Senator, 
DeCarlo,  and  Torrillo  (Ref.  27) 
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Comparison  of  Present  Results  with  Experimental  Data  of  Wilson 
and  Fisher  (Ref.  28)  and  Wilson  (Ref.  29) 
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TABLE  I 

SUMMARY  OF  EXPERIMENTAL  DATA  USED  IN  COMPARISONS 
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P0>®J  psia 

To  ,  °R 

L»,  in. 
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TABLE  II 

SUMMARY  OF  INVISCID  EDGE  CONDITIONS  USED  IN  COMPARISONS 
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APPENDIX  III 

IMPLICIT  FINITE-DIFFERENCE  SOLUTION  OF 
GOVERNING  BOUNDARY-LAYER  EQUATIONS 
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For  the  special  case  of  a  sharp  cone  in  a  hypersonic  flow,  the  invis- 
cid  flow  field  is  conical  in  character  so  that  all  the  inviscid  flow  param¬ 
eters  (p,  Te,  pe,  Ue)  remain  constant  along  rays  from  the  apex  of  the 
cone.  Under  this  restriction  the  governing  equations  (43)  and  (44)  in  the 
transformed  plane  reduce  to  the  following  simpler  set 
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which  apply  to  any  boundary -layer  flow  having  constant  inviscid  flow 
parameters.  The  boundary  conditions  remain  unchanged  from  Eqs.  (52) 
through  (56)  which  are  repeated  below  for  completeness. 


MOMENTUM 


f(5,  n  =  0)  =  0 

(III-3) 

f';(5,  n  =  0)  =  0 

(III-4) 

lim  f'(5»n)  *  1 

n  » 

(IH-5) 
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ENERGY 


e(e«  n  =  0)  =  ew  U) 


(in-6) 


lin  9(c,n)  =  1 
n  -*■  • 


(III -7) 


Following  the  approach  by  Davis  (Ref.  43)  and  Blottner  (Ref.  44)  the 
momentum  and  energy  equations  (Eqs.  (Ill- 1)  and  (III-2))  are  rewritten 
in  "standard"  form  for  a  parabolic  partial  differential  equation  as 

a2W  ^  3W  .  „  3W  _  n 

3n2  +  “1  3n  “2”  +  “3  +  «4  3?  "  (III-8) 


where  W  =  f'for  the  momentum  equation  and  W  =  9  for  the  energy  equa¬ 
tion.  Using  Eqs.  (HI-1)  and  (111-2),  the  coefficients  cq  through  04  are 
found  to  be,  in  linearized  form: 


MOMENTUM 


a  a 


-  in 


-  +  f  + 


a* 


^  as 


“2 

“3 


=  0 

=  0 


-  25  f' 
"  2* 


ENERGY 


a 

a** 

,  „  3  f 

an 

Pr  ' 

tf  +  25  at 

Pr 

a2  =  0 


(III- 9) 
(III- 10) 
(HI- 11) 
(IH-12) 

(III- 13) 

(III- 14) 
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(m-i5) 


(in- 16) 


The  r)  derivatives  in  Eq.  (III-8)  are  replaced  with  finite -difference 
quotients  which  allow  variable  grid  spacing  in  the  rj  direction  so  as  to 
concentrate  grid  points  in  the  region  near  the  body  surface  where  the 
dependent  variables  change  most  rapidly.  The  derivative  in  the  ?  di¬ 
rection  in  Eq.  (III-8)  is  handled  in  the  usual  manner  as  a  two-point 
backward  difference  between  points  (m+l,n)  and  (m,n)  while  all  r\ 
derivatives  are  evaluated  at  point  (m+1,  n)  according  to  the  grid  mesh 
shown  below. 


Station 


Boundary-  m  m+1 


The  solution  is  assumed  to  be  known  at  point  (m,  n)  and  unknown  at 
point  (m+l,n)  so  that  the  finite -difference  scheme  to  be  constructed  will 
be  implicit  in  nature.  The  finite -difference  replacements  for  the  deriva¬ 
tives  are  as  follows: 


a2wn 


drj  2 


2[  W  +  KW  ,  -  <1+K)W  1  ^ 
=  L  n+1  n-1  nJm+l 


■m+1,  n 


TX 


(III- 17) 
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IS]  =  t  Wn+1  -  K2Wn.!  -  (l-K2)Wn5  m+1 
^  Jm+l,n 


(in- is) 


3>j]  _  wm+1,n 
35  ~ 


tyn+Un  ~  Wf,i,n 


m+1  ,n 


(III- 19) 


where 


D1  =  K+l  “  %)  +  '<2(nn  -  nn_1 ) 


(in- 20) 


d2  s  (nn+l  "  nn)2  +  KK  "  nn-l)' 


(in-21) 


K  =  -  —  (constant) 

nn  "  nn-1 


(in-22) 


=  n2  "  nl  (constant) 


(in-23) 


The  finite -difference  form  of  Eq.  (UI-8)  becomes,  upon  substitution  of 
Eqs.  (in-17),  (m-18),  and  (111-19), 

^n  tyn+1  ,n+l  +  ^m+1  ,n  +  ^n  Hn+1  ,n-l  =  (ni-24) 


where 


.  .  2  al 
An  ‘  +  Ef 


(m-25) 


-2(1+K)  «i (1 -K2 ) 


+  a2  +  A? 


(HI- 2  6) 


2_K  ^2°1 

^n  ^2  t'-j 


(HI- 2  7) 


Rn  =  "  a3  + 


a4  n 


(in-28) 


In  order  for  Eq.  (DI-24)  to  be  linear,  the  coefficients  ■^n*  ®n*  ^n» 
and  Rn  must  be  treated  as  known  quantities  at  point  n;  more  will  foHow 
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on  this  subject  later.  The  important  point  is  that  Eq.  (IH-24)  represents 
a  set  of  simultaneous  linear  algebraic  equations  under  this  restriction. 

Since  the  simultaneous  linear  algebraic  equations  resulting  from 
Eq.  (III-24)  are  of  a  special  form  (tridiagonal),  an  efficient  method  of 
solution  on  a  digital  computer  is  available  from  Richtmyer  and  Morton 
(Ref.  45,  pp.  198  to  201  and  274  to  282).  For  this  procedure  the  bound¬ 
ary  conditions  at  the  wall  (n  =  1)  and  the  outer  edge  (n  =  N)  must  have 
specified  values  i  and  Wm+ij  jj.  Because  of  the  special  form  of 

Eq.  (III-24),  the  relation 


wm+l ,n  =  En  wm+l  ,n+l  + 

eir» 

2  -  n  -  N-1 

(HI- 2  9) 

exists  where 

c  ,± 
e2 

(HI-30) 

p2  -  c2  wr*1 ,1 

(IH-31) 

e2  "  E>2 

r  _  _An  N 

tn  -  r  - +  r  t — r 

Ln  n-1 

(III- 3  2) 

„  _  Rn  "  cn  en-l  ! 

n  Bn  +  cn  En-1 

3  -  n  -  N-1 

(III- 3  3) 

The  quantities  En  and  en  are  computed  from  Eqs.  (HI-30)  through 
(HI-33)  starting  with  n  =  2  and  progressing  to  n  =  N-1.  The  solution 
wm+l,n  is  then  obtained  by  evaluating  Eq.  (IH-29)  from  n  =  N  -  1  to 

n  =  2.  Knowing  the  distribution  of  f '  and  6  across  the  boundary  layer 
from  the  above  procedure,  the  transformed  stream  function  f  is  eval¬ 
uated  from 


f(?,n)  = 


f'  (S,n)dn 


(III-34) 
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where  the  integral  is  numerically  integrated  using  the  well-known  trape¬ 
zoidal  rule,  viz 

f(£.n  =  nn)  -  ff'U.rt  8  nj  )  +  f'U.n  *  ni.1)joi  (ni-35) 
1=2 


with 


Di 


"1 


-  n1-l 
2 


(HI-36) 


In  a  similar  manner,  the  inversion  from  the  transformed  (5,rj)  plane  to 
the  physical  (x, y)  plane  is,  from  Eq.  (37), 


where 


y  = 


M. 


feVb0 


/ 


pe  dn 

p" 


(in-37) 


(HI-38) 


because  of  the  constancy  of  static  pressure  across  the  boundary  layer. 
Again,  using  the  trapezoidal  rule  method  of  numerical  integration 
yields 


y 


n 


peW 


0(?*n=  n-j )  +  e(C*  n=n1-_-] 


(III- 3  9) 


with  Dj  given  by  Eq.  (HI-36).  The  relationship  between  I  and  x  for 
sharp  cone  flows  is  found  by  integration  of  Eq.  (36)  with  rfcJ  =  x  sin  6V 
from  geometry  which  results  in 


5  = 


1 

3 


pevel)e  sin2  «v  x: 


(IH-40) 


For  the  present  work  the  A?  increment  of  integration  is  chosen  based  on 
the  division  Ax  =  0.  01  L,  i.  e. ,  100  stations  spaced  equally  along  the 
physical  body. 


The  mathematical  basis  of  the  above  tridiagonal  matrix  procedure 
applied  to  the  solution  of  boundary -layer  problems  is  attributable  to 
Flugge-Lotz  and  HLottner  (Ref.  46).  The  present  application  differs 
from  their  original  work  in  one  important  aspect:  the  linearized  dif¬ 
ference  equations  herein  are  uncoupled  and  solved  separately.  In 
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Flugge-Lotz  and  Blottner's  approach,  the  difference  equations  remain 
coupled  and  require  additional  machine  storage  and  manipulations  for 
solution.  With  the  present  uncoupled  approach,  the  difference  equations 
are  iterated  to  convergence  at  each  station  along  the  body  so  that  one 
must  pay  the  price  of  iteration;  this  procedure  is  how  the  linearizing 
coefficients  aj  through  <24  are  evaluated  at  each  station  using  the  results 
of  the'  previous  iteration.  Hence  the  final  solution  obtained  at  each 
station  is  exact  in  the  sense  that  it  represents  a  converged  iterated  solu¬ 
tion  to  the  governing  nonlinear  partial  differential  equations  written  in 
finite -difference  form.  Typically  about  three  iterations  per  station  are 
required  for  fully  turbulent  flow;  only  one  to  two  iterations  per  station 
are  necessary  for  laminar  flow. 

The  variable  grid  mesh  used  in  the  present  work  is  taken  from 
Smith  and  Cebeci  (Ref.  2).  The  various  constants  used  herein  are  as 
follows: 

N  =  85 
K  =  1.063 
Q  =  0.010 

Experience  with  varying  these  constants  and  observing  their  influence  on 
the  resultant  numerical  solution  has  indicated  that  the  above  choices  are 
adequate  under  the  present  flow  conditions.  These  values  may  not  be 
satisfactory,  however,  for  other  body  geometries  and  flow  conditions  so 
that  the  influence  of  the  variable  grid  mesh  constants  should  be  ascer¬ 
tained  for  each  new  investigation.  Provisions  are  included  in  the  pres¬ 
ent  program  for  addition  of  points  (if  necessary)  to  the  solution  as  the 
body  is  traversed.  This  may  be  required  in  order  to  prevent  numerical 
suppression  of  boundary -layer  growth. 

At  the  apex  of  a  sharp  cone  where  §  =  0,  the  governing  boundary - 
layer  equations  (Eqs.  (Ill—  1)  and  (III— 2))  reduce  to  the  following  ordinary, 
nonlinear,  differential  equations 


MOMENTUM 


Z*  f 


f" 


=  0 


(III-41) 
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ENERGY 


+  (y-1)!^2  A*  (f")2  -  0 


(in- 42) 


with  the  boundary  conditions  of  Eqs.  (Ill- 3)  through  (III- 7).  In  order  to 
obtain  starting  profiles  for  the  finite -difference  scheme  to  march  down¬ 
stream,  Eqs.  (IH-41)  and  (HI-42)  are  solved  by  using  the  tridiagonal 
matrix  procedure  described  previously  in  the  foUowing  manner.  At 
the  apex  of  the  cone  §  =  0  so  that  04  =  0  in  both  the  momentum  and 
energy  equations;  in  addition,  the  term  containing  §  in  <*1  of  both  the 
momentum  and  energy  equation  vanishes*.  Initial  guesses  for  f'  and  0 
are  input  to  the  program  as 

eI 

where  the  subscript  I  denotes  the  initial  approximation.  The  equations, 
are  then  iterated  to  convergence  in  the  same  manner  described  previ-J  , 
ously;  an  averaging  scheme  is  used  to  speed  convergence.  Typically; 
about  20  to  25  iterations  are  required  to  generate  a  converged  initial 
solution. 


■  1  -  exp  (-n) 


:  -  ew  +  0  -  ew)fi 


(HI-43) 
(IH- 44) 


By  use  of  the  above  procedures,  the  numerical  solution  of  any  two- 
point  boundary -value  problem  governed  by  either  linear  or  nonlinear 
ordinary  differential  equations  as  weH  as  sets  of  coupled  parabolic 
partial  differential  equations,  either  linear  or  nonlinear,  is  reduced 
to  subroutine  status  on  a  digital  computer  in  that  only  the  coefficients 
a 1  through  a 4  must  be  defined  in  conjunction  with  the  required  boundary 
conditions  for  each  new  problem.  Much  application  of  these  procedures 
has  been  found  in  the  von  Karman  Gas  Dynamics  Facility  of  the  Arnold 
Engineering  Development  Center  on  many  different  problems  involving 
both  nonlinear,  ordinary,  differential  equations  as  well  as  parabolic, 
partial,  differential  equations.  The  thin  viscous  shock  layer  analysis 
by  Adams  (Ref.  47)  is  a  good  example.  In  this  work  eight  simultaneous, 
nonlinear,  ordinary,  differential  equations  governing  momentum, 
energy,  and  species  conservation  were  solved  using  the  iterative  tri¬ 
diagonal  matrix  method  presented  above.  The  problem  was  made  even 
more  difficult  in  that. chemical  nonequilibrium  effects  were  included  in 
the  analysis.  Based  on  experience  with  analyses  of  this  type,  the 
author  highly  recommends  use  of  the  iterative  tridiagonal  matrix 
approach  where  applicable. 
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Analysis  of  laminar,  transitional,  and  turbulent  boundary  layers  on 
a  sharp  cone  at  zero  angle  of  attack  under  hypersonic  perfect  gas  condi¬ 
tions  is  presented.  The  governing  boundary-layer  equations  were 
numerically  integrated  on  a  digital  computer  using  a  marching,  iterative, 
implicit,  finite-difference  method.  The  turbulent  boundary  layer  is 
analyzed  using  a  two-layer  eddy  viscosity  model  with  a  constant  turbu¬ 
lent  Prandtl  number,  and  the  transition  region  is  treated  through  an 
eddy  viscosity-intermittency  factor  approach.  Comparison  with  experi¬ 
mental  data  reveals  that  the  present  approach  yields  an  accurate  de¬ 
scription  of  laminar,  transitional,  and  fully  turbulent  heat  transfer  on 
sharp  cones  in  the  Mach  number  range  from  5  to  10  under  cold  wall  con¬ 
ditions.  Verification  of  the  detailed  calculation  of  transitional  and 
turbulent  boundary-layer  structure  under  hypersonic  conditions  must 
await  further  experimental  investigations.  It  appears  at  the  present 
time  that  the  eddy  viscosity  approach  to  calculation  of  transitional  ( 

and  turbulent  boundary  layers  may  indeed  be  applicable  and  accurate  ‘ 

under  hypersonic  conditions.  ' 
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